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TECHNICAL  SUMMARY 


The  major  thrust  of  the  research  accomplished  under  this  grant  ia 
toward  algorithms  that  extend  the  size  of  problems  we  can  effectively 
solve*  although  there  have  been  other  results  along  the  way*  -Most  of 
this  technical  material  has  been  previously  reported*  but  there  has  been 
substantial  progress  since  the  last  report  and  we  will  concentrate  on 
that  here. 

This  project  began  with  the  main  goal  of  finding  a  variable-metric 
technique  for  large  sparse  nonlinear  optimization  problems  by  the  novel 
approach  of  defining  the  sparse  derivative  matrix  approximations  in 
terms  of  their  sparse  factorizations  which  would  be  updated  at  each 
iteration.  We  are  preparing  a  paper  dealing  with  our  findings  for  the 
Cholesky  factorization.  Our  idea  seems  to  be  sound*  but  more  work  is 
needed  to  make  the  work  practical.  On  the  other  hand*  the  LU  form* 
introduced  in  the  joint  paper  with  Marwil  for  nonlinear  simultaneous 
equations*  is  a  success.  An  interesting  extension  of  that  work  by 
researchers  in  Brazil  is  forthcoming.  Most  people  in  optimization  have 
now  given  up  on  the  more  conventional  approach  to  sparse  variable-metric 
methods  because  of  the  scale  dependence  of  the  sparsity  property,  so  we 
expect  our  factor-updating  approach  to  be  considered  more  in  the  future. 

Franklin  Luk  has  been  working  with  Professor  Richard  Brent  of  the 
Australian  National  University  on  algorithms  for  systolic  array 
machines.  This  collaboration  has  so  far  produced  four  joint  papers* 
copies  of  which  are  enclosed.  These  papers  continue  work  previously 
described  in- progress  reports  on  the  exploitation  of  this  novel  machine 
architecture  in  parallel  numerical  linear  algebra  computations. 

Since  the  last  progress  report,  John  Dennis  and  Homer  Walker  have 
obtained  and  written  up  the  results  on  local  improvement  theorems  which 
are  enclosed  with  this  report.  It  is  always  difficult  to  tell  how 
results  will  be  received  by  the  numerical  community*  but  these  results 
seem  to  be  a  great  help  in  understanding  the  effects  of  modeling  and 
computational  errors  in  determining  the  attainable  accuracy  and  speed  of 
convergence  to  a  solution  of  that  accuracy  for  a  system  of  nonlinear 
equations.  The  final  goals  of  this  research  are  to  furnish  a  front-end 
to  our  research  reported  in  a  previous  progress  report  on  exploiting 
*almost-sparsity'  in  nonlinear  problems,  especially  discretizations  of 
nonlinear  pde's*  and  more  fundamentally,  we  hope  eventually  to  modify 
library  subroutines  like  MINPACK  to  supply  the  user-supplied  function 
subroutine  with  an  accuracy  requirement  with  each  function  evaluation 
request.  This  has  the  potential  for  enormous  computational  savings  in 
inverse  problems  where  the  function  evaluation  is  the  result  of  a  numer¬ 
ical  simulation,  since  it  would  furnish  a  rational  scheme  for  a  less 
accurate  simulation  early  in  the  parameter-identification  process. 
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